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Preface 


In  1980,  the  C-141B  modification  program  began  in  earnest. 
Suddenly,  C-141  flight  crews  had  to  develop  expertise  in  an  area 
in  which  they  had  virtually  no  experience  -  air  refueling.  The 
difficulties  which  many,  including  myself,  encountered,  prompted 
this  research.  Above  all,  this  study’ s  goal  is  simply  to  make 
things  easier  for  the  crewmember,  who  is  given  a  variety  of 
missions  and  limited  training  resources.  Hopefully,  this  paper 
will  cause  others  to  look  more  closely  at  the  air  refueling 
problem  in  order  to  develop  simpler,  more  efficient  methods. 

I  am  indebted  to  many  people  for  their  help  in  this 
effort.  I  am  especially  grateful  to  Lt  Col  Edwards  for  all  of 
his  advice,  guidance,  and  patience.  I  would  also  like  to 
thank  all  of  those  who  assisted  in  the  research,  particularly 
Lt  Ham  and  Greg  Santaromita  at  Wright-Patterson,  Maurice  Johns 
with  Delco,  and  Bob  Harper  at  Collins.  And,  I  am  grateful  that 
the  typist,  Linda  Kankey,  can  meet  a  deadline  better  than  I  can. 


James  T.  Rivard 
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Abstract 

XsH.  Estimation  and  control  algorithms  were  developed  for 
use  by  a  tanker  aircraft  conducting  an  air  refueling  rendezvous. 

A  stochostic  model  of  a  typical  rendezvous  was  developed  first. 
Then  an  entended  Kalman  filter  which  uses  air-to-air  TACAN 
distance  measurements  was  designed.  Also,  algorithms  were 
derived  for  computing  the  tanker's  appropriate  offset,  turn 
point,  and  closed  loop  bank  angle  commands  during  the  final 
turn  of  a  point-parallel  rendezvous - 

In  Monte  Carlo  simulations,  a  3  state  extended  Kalman 
filter  provided  a  stable,  though  limited,  means  of  estimating 
the  receiver's  position  and  velocity  throughout  the  rendezvous. 
Also,  the  control  algorithms  exhibited  two  advantages  over  the 
present  rendezvous  method:  turn  point  solutions  could  be  computed 
for  other  than  nominal  offsets  and  headings,  and  bank  angle 
commands  during  the  turn  reduced  position  errors  at  rendezvous 
completion.  When  compared  to  the  present  method,  root  mean 
square,  cross-track  distance  errors  at  rollout  were  reduced  as 
much  as  75$.  and  along-track  errors  were  reduced  up  to  47*. 
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DIGITAL  ESTIMATION  AND  CONTROL 
FOR  AIR  REFUELING  RENDEZVOUS 


I  Introduction 

Certainly,  in-flight  refueling  is  a  vital  element  of  Air 
Force  strategy  and  tactics.  It  has  long  been  a  critical  part  of 
bomber  and  fighter  operations.  Recently  it  has  assumed  an 
increasingly  important  role  in  airlift  planning  as  well,  parti¬ 
cularly  with  the  advent  of  the  Rapid  Deployment  Force.  The  Air 
Force's  large  fleet  of  KC-135  and  KC-10  aircraft  provide  the 
United  States  with  a  unique  military  capability.  No  other  nation 
in  the  world  can  project,  in  a  matter  of  hours,  such  a  wide  range 
of  military  responses  to  a  perceived  threat  anywhere  in  the  world. 

Of  course  the  great  advantage  of  in-flight  refueling  is 
that  it  eliminates  the  need  for  enroute  airfields.  Many  naviga¬ 
tional  systems  have  been  designed  to  get  an  aircraft  safely  into 
an  airfield.  However,  there  are  no  systems  in  operational  use 
designed  specifically  for  conducting  air  refueling  rendezvous. 
Instead,  a  wide  variety  of  navigation  and  radar  systems  designed 
primarily  for  other  purposes  are  adapted  to  the  problem.  The  two 
aircraft  usually  get  together, but  the  initial  rendezvous  attempt 
is  often  very  inaccurate.  That  frequently  causes  excessive 
maneuvering  and  delays  in  completing  the  join-up.  Those  undesir¬ 
able  effects  take  on  additional  significance  for  large  aircraft 
requiring  extended  refueling  time,  aircraft  formations,  and 
refueling  operations  at  night  or  in  instrument  weather  conditions. 
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Hardware  technology  exists  which  should  meet  any  practical 
rendezvous  specification.  The  first  step  in  upgrading  our  air 
refueling  systems,  though,  is  to  investigate  what  improvements 
can  be  made  with  the  existing  equipment. 

The  KC-135  fleet  was  recently  equipped  with  an  inertial 
navigation  system  (INS).  An  integral  part  of  the  INS  is  the 
navigation  computer.  This  computer  uses  information  from  the 
inertial  instruments  and  the  central  air  data  computer  (CADC)  to 
compute  aircraft  position  and  velocity,  and  wind  information. 

It  also  calculates  steering  commands  for  the  autopilot  and  flight 
director  system.  The  KC-135  and  all  Air  Force  refuelable  aircraft 
are  also  equipped  with  Air-to-Air  TACAN  (Tactical  Air  Navigation). 
The  system  provides  a  measurement  of  straight  line  distance  to 
other  TACAN  equipped  aircraft. 

A  possible  system  improvement  would  be  for  the  tanker  to 
use  TACAN  measurements  as  inputs  to  an  estimation  algorithm  which 
has  outputs  of  the  other  aircraft's  position  and  velocity.  These 
estimates  would  then  be  used  to  compute  steering  commands  for  the 
tanker  during  the  final  part  of  the  rendezvous.  If  all  of  those 
estimation  and  guidance  computations  can  be  handled  by  the  present 
onboard  navigation  computer,  the  system  could  be  implemented 
entirely  on  current  hardware. 

The  purposes  of  this  study  are  to  develop  an  estimation 
and  control  method  for  solving  the  rendezvous  problem  and  to  see 
what  improvements  in  accuracy  such  a  system  might  achieve. 


Along  with  the  existing  hardware  constraint,  operational  feasi¬ 
bility  is  a  prime  consideration  throughout.  The  standard 
rendezvous  profile  is  altered  as  little  as  possible  so  that 
present  methods  can  always  be  used  as  a  back-up. 
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The  most  common  type  of  rendezvous  is  called  the  point 
parallel  rendezvous.  In  this  type,  the  tanker  waits  for  the 
receiver,  the  aircraft  to  be  refueled,  at  a  planned  point  along 
the  receiver's  route  of  flight.  The  rendezvous  is  accomplished 
and  the  refueling  completed  on  a  refueling  track  in  the  receiver's 
general  direction  of  flight.  This  type  of  rendezvous  will  be 
discussed  here. 

Rendezvous  Problem 

A  typical  point  parallel  rendezvous  is  depicted  in  Figure 
1 .  The  tanker  waits  at  refueling  altitude  in  a  holding  pattern 
at  the  control  point  (CP).  When  the  receiver  reports  that  he  has 
passed  the  initial  point  (IP),  the  tanker  departs  the  holding 
pattern  and  flies  towards  the  receiver.  Meanwhile,  the  receiver 
proceeds  along  the  refueling  track  at  constant  airspeed  and 
descends  to  an  altitude  1000  feet  below  the  tanker.  The  tanker 
attempts  to  establish  the  correct  offset  and  turn  at  the  right 
time  so  as  to  roll  out  of  a  30  degree  bank  turn  two  nautical 
miles  directly  in  front  of  the  receiver.  The  receiver  then 
closes  to  one  mile  from  which  point  the  rendezvous  is  completed 
visually. 

The  rendezvous  problems  for  the  tanker,  then,  are  to 
compute  and  fly  the  appropriate  offset,  turn  at  the  correct  time, 
and  maintain  proper  bank  angle  throughout  the  turn. 
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Figure  1,  Point  Parallel  Rendezvous 
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Present  Method 


Currently,  a  chart  is  used  to  determine  the  correct  offset 
and  turn  point.  Given  values  for  receiver  true  airspeed  and 
drift  angle,  the  chart  gives  the  appropriate  offset  and  distance 
from  the  receiver  to  initiate  a  30  degree  bank  turn.  Here, drift 
angle,  the  number  of  degrees  the  aircraft  must  crab  into  the  wind 
to  maintain  a  course  parallel  to  the  refueling  track,  is  used  as 
a  measure  of  the  wind  component  perpindicualr  to  the  track. 

Since  the  tanker  and  receiver  navigation  systems  do  not  agree 
exactly,  they  attempt  to  establish  the  desired  offset  by  observing 
the  other  aircraft's  position  on  radar.  When  the  tanker's  TACAN 
and/or  radar  indicates  the  receiver  is  at  the  calculated  turn 
range,  the  tanker  rolls  into  a  30  degree  bank  turn  which  he  holds 
until  he  has  turned  to  the  refueling  heading. 

The  present  method  has  several  inaccuracies.  The  offset 
and  turn  range  chart  only  gives  distances  to  the  nearest  mile; 
fluctuations  in  wind  or  airspeed  are  not  compensated  for;  equipped 
with  a  radar  designed  primarily  for  weather,  it  is  difficult  for 
the  KC-135  to  determine  accurately  another  aircraft's  position; 
and,  finally,  the  system  is  entirely  open  loop;  a  constant  30 
degrees  of  bank  is  held  throughout  the  turn.  Even  under  the 
best  visual  conditions  it  is  difficult  for  either  the  tanker 
or  receiver  to  recognize  any  need  for  corrections  until  the  turn 
is  nearly  complete. 
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Proposed  Method 


0> 


The  proposed  method  would  be  basically  the  same  proce¬ 
dure  but  with  two  improvements.  Rather  than  using  a  chart,  the 
correct  offset  and  turn  point  would  be  calculated  and  periodically 
updated  by  the  navigation  computer.  Also,  an  estimation  algorithm 
using  Air-to-Air  TACAN  distance  measurements  would  be  used  to 
estimate  continuously  the  receiver's  position  and  velocity. 

That  estimate  would  be  used  as  input  to  a  closed  loop  digital 
control  algorithm  before  and  during  the  turn. 

A  rendezvous  conducted  with  such  a  system  might  proceed 
as  follows: 

1.  Rendezvous  data  is  loaded  into  tanker's  INS 
anytime  prior  to  rendezvous.  Data  includes 
control  point  coordinates,  refueling  track 
heading  and  altitude,  and  receiver  altitude 
and  airspeed. 

2.  Using  these  input  data  and  information  from  the 
tanker' s  INS  and  central  air  data  computer, 

the  navigation  computer  calculates  and  displays 
the  appropriate  offset. 

3.  Tanker  and  receiver  attempt  to  establish  correct 
offset  using  radar. 

4.  Estimator  is  initialized  based  on  planned 
receiver  track,  altitude,  and  airspeed,  and 
computed  wind  velocity.  A  manual  radar  update 
might  also  be  used. 
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5.  Estimation  algorithm  propagates  receiver  state 
estimate  and  uses  TACAN  measurements  for  updates 

6.  Computer  calculates  and  displays  time  until  turn 
and  commands  a  turn  at  updated  turn  point. 

7.  Digital  controller  provides  bank  angle  commands 
to  the  autopilot/flight  director  to  reduce  error 
at  rollout. 

The  first  step  in  designing  such  a  system  is  to  create 
a  mathematical  model  of  the  rendezvous  procedure  just  described. 
That  model  will  later  be  used  as  the  starting  point  for  the 
estimator  design,  and  as  a  reference  or  truth  model  for  computer 


simulations. 


Ill  Truth  Model 


Research  was  conducted  into  present  air  refueling  proce¬ 
dures,  the  KC-135  inertial  navigation  system,  TACAN  systems, 
and  wind  models.  The  result,  as  developed  in  this  section,  is 
a  stochastic  state  differential  equation  and  measurement  equation 
which  hopefully  describe  the  actual  dynamics,  uncertainties,  and 
errors  reasonably  well.  For  different  types  of  receiver  aircraft, 
the  form  of  the  state  equations  will  be  the  same.  However,  there 
will  be  slight  variations  in  certain  parameters  and  error  statis¬ 
tics  due  to  different  rendezvous  airspeeds  and  navigation  equip¬ 
ment.  In  the  truth  model  development  and  the  simulations,  a 
standard  rendezvous  with  a  C-141  is  modelled.  This  study,  while 
limited  to  one  aircraft  type,  should  indicate  both  the  practical- 
ity  and  the  potential  benefit  of  the  proposed  system. 

Coordinate  System 

The  coordinate  system  for  the  truth  model  is  depicted 
in  Figure  2.  It  is  a  simple,  two-dimensional,  rectangular 
system  parallel  to  the  tanker's  local  horizon  plane.  The  origin 
is  located  at  the  control  point  (CP)  at  refueling  altitude.  The 
plus  x  axis  points  towards  the  initial  point  (IP),  and  the  plus 
y  axis  is  oriented  as  shown.  Aircraft  headings,  T,  are  measured 
counter-clockwise  from  the  plus  x  axis. 


Figure  2.  Coordinate  System 


The  exact  reference  frame  location  and  orientation  are  defined 
by  the  INS,  so  the  frame  varies  slightly  from  a  true  ground 
fixed  system  due  to  INS  errors.  The  advantage  of  defining  the 
system  by  the  INS  is  that  INS  tanker  position  readouts  are 
deterministic  rather  than  random  variables. 

Rectangular  coordinates  are  chosen  to  make  the  filter 
state  propagation  equations  linear.  Also,  this  particular  axes 
orientation  simplifies  turn  point  calculations:  y  position 
coordinates  determine  offset,  and  x  position  and  velocity  com¬ 
ponents  determine  the  time  to  turn. 

Receiver  Dynamics 

For  an  aircraft,  the  velocity  in  a  ground  fixed  frame  is 
found  from 


*  =  — a/w  T 
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where  v  is  velocity  relative  to  the  ground,  v  /  is  velocity  of 

the  aircraft  relative  to  the  wind  or  air  mass,  and  v  is  the 

— w 

velocity  of  the  wind. 

The  receiver  attempts  to  fly  a  ground  track  along  the 
x  axis.  In  order  to  do  so,  he  must  apply  a  drift  correction,  6, 
into  the  wind  such  that  the  y  component  of  velocity,  v  ,  is  zero. 

y 

Since  he  is  heading  in  the  -x  direction,  his  heading  will  then 
be  ¥=tt  +  6. 

The  receiver  pilot  will  also  attempt  to  maintain  a 
preplanned  indicated  airspeed.  The  tanker  can  convert  this 
indicated  airspeed  to  an  equivalent  true  airspeed,  v^,  based  on 
rendezvous  altitude  and  measured  air  temperature.  The  components 
of  are  then 

va/w  x  "  VR  cos  (w  +  41 


va/„  „  =  Vn  sin  (tt  +  <$) 


vSl/ w  y  R 

or,  equivalently. 


va/ w  x  *  ”VR  005  S 


v„ /  =  -vD  sin  5 

a/w  y  R 

If  the  receiver  maintains  v  =  0,  the  y  component  of  Eq ( 1 ) 


becomes 


0  =  -vR  sin  6  +  vwy 

which  can  be  used  to  find  <5. 
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The  navigation  computer  continuously  calculates  the  wind  vector, 
vw»  by  subtracting  va/w“  determined  from  measured  true  airspeed 
and  aircraft  heading  -  from  the  INS  computed  velocity,  v.  The 
planned  refueling  track  heading  can  then  be  used  to  calculate 

v„x  and  vwy’ 

Now  the  velocity  relations  for  the  receiver  can  be 
written.  The  receiver  attempts  to  fly  such  that 

6  =  sin  "  [  vwyl  (: 


v  =o 

y 


vx  *  -VR  cos  6  +  v„x 

Of  course  his  flight  will  not  exactly  satisfy  those  equations. 
Instead,  the  actual  velocity  relations  will  be 


v  =  v 
y  ye 


v  =  - vn  cos  6  +  v  +  v 
x  R  wx  xe 

where  v  and  v  are  the  velocity  error  components, 
xe  ye 

Cross-Track  Error 

As  mentioned  previously,  the  modelling  parameters  in 
this  study  are  based  on  a  standard  rendezvous  of  a  KC-135  and  a 
C-141.  Both  aircraft  are  equipped  with  a  Carousel  IV  inertial 
navigation  system  made  by  Delco.  The  C-141' s  autopilot  can  very 
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accurately  maintain  the  course  computed  by  its  INS.  The  INS, 
though,  drifts  slightly  during  flight  and  is  the  source  of  some 
cross-track  error. 

Data  from  several  hundred  military  and  commercial  flights 
indicated  a  95$  confidence  figure  of  1.65  NM/hr  for  the  Carousel 
IV  INS.  According  to  Delco,  a  95%  cross-track  confidence  value 
can  then  be  estimated  as  66$  x  1.65  NM/hr  =  1.1  NM/hr  (Ref  1:  6-1). 
Using  this  as  a  2 a  value,  the  INS  has  ala  cross- track  error  rate 
of  0.55  NM/hr. 

Position  errors  in  an  INS  have  a  variety  of  sources  and 
the  associated  error  models  can  be  very  elaborate  (Ref  2) . 

However,  for  the  short  time  interval  of  a  rendezvous  -  less  than 
10  minutes,  the  cross-track  error  might  be  reasonably  modelled 
as  the  sum  of  a  constant  rate  component  and  a  small  random  Walk 
component.  Because  the  rendezvous  duration  is  brief  compared 
to  the  84  minute  Schuler  period,  and  because  this  study  is  for 
feasibility  and  not  detailed  operational  design,  second-order, 
coupled  error  models  are  not  included. 


The  random  constant  error  rate,  v  ,  models  the  long  term 

y  c 

position  error  source  which  flight  tests  have  shown  to  be  the 
dominant  source  of  velocity  error  for  the  Carousel  IV  (Ref  1:  6-6). 
VyC  is  assumed  to  be  Gaussian  and  zero  mean,  and  to  have  a  standard 
deviation  of  0.55  NM/hr.  It  is  also  assumed  to  be  stationary: 


the  error  rate  over  a  short  period  has  the  same  statistics  as 
the  error  rate  over  the  whole  flight.  But  that  is  only  the 
error  rate  due  to  the  receiver’s  INS.  Remember  that  the  receiver's 


position  and  velocity  are  estimated  in  the  reference  frame  of 
the  tanker's  INS.  Any  error  in  the  tanker's  INS  is  then  perceived 
as  an  error  in  the  receiver's  state.  So  the  total,  apparent, 
constant  drift  rate  is 


v  =  v  _ 
yc  yc  Receiver 


yc  Tanker 


Since  each  INS  operates  independently,  and  the  two  aircraft  may 
fly  entirely  different  routes  to  the  rendezvous,  their  random 
drift  rates  are  assumed  to  be  independent.  Then  v^c  is  zero 
mean,  Gaussian,  and 


a  =  a  r.  • 
vyc  vyc  Receiver 


+  a 


vyc  Tanker 


avyC  is  thus  calculated  to  be  0.778  NM/hr. 

A  random  walk  component  of  position  error,  y^,  is  added 

to  account  for  small  deviations  from  the  simplistic,  constant 

error  rate  model.  This  component  is  modelled  by 

y  (t)  =  w  (t) 

'wv  y 

where  w  (t)  is  a  zero  mean,  continuous,  white  noise,  and 

y 

E{Wy(t)  Wy  (t  +  T)}  =  <ly  5  (  T ) 

E{»}  is  the  expected  value  operator,  6(*)  is  the  Dirac  delta 
function,  and  q  is  called  the  strength  of  the  white  noise. 

v 

This  type  of  model  results  in  a  zero  mean  process  with 


E  (yw  (t) }  =  qy(t-tQ) 

where  t  is  the  initial  time  of  the  simulation, 
o 


14 


Choosing  an  appropriate  value  for  q  was  somewhat 

y 

arbitrary.  Schuler  oscillations  contribute  a  0.2  NM/hr  CEP 
velocity  error  component  (Ref  1:  6-6).  Over  the  short  10  minute 
rendezvous  period,  a  0.2  NM/hr  error  would  cause  a  202  ft  position 
error.  To  account  for  both  aircraft’s  errors,  INS  errors  for 
which  no  error  data  could  be  found,  and  autopilot  errors,  a  1o 


value  for  y  at  the  end  of  a  10  minute  period  was  chosen  to  be 
500  ft.  The  idea  here  was  to  choose  a  conservative  (high)  value 
for  the  magnitude  of  the  errors  not  modelled  by  vyc.  Since  the 
proposed  method  would  have  to  be  suitable  for  several  receiver 
aircraft  types,  a  specific  C— 1 41  INS  error  model  is  inappropriate. 
The  noise  strength  is  then  found  to  be 


qy  = 


(500  ft) 
(10  min) 


=  25000  ft  /min 


5  2 

=  1.13  x  10"  NM  /sec 

The  model  for  cross-track  position  and  velocity  is  now 
written  as 

y(t)  =  vye(t) 

vye(t)  *  V  *  “y(t) 

E{vyc}  =  0,  avyc  =  2.161  x  10"  NM/sec 

5  2 

E{ w  (t)}  =  0,  q  =  1.13  x  10"  NM  /sec 

y  y 


Now,  only  the  initial  condition  statistics  for  y(t)  remain 
to  be  specified.  In  the  simulations,  two  different  situations 
are  modelled.  In  one  case,  the  rendezvous  is  conducted  without 
the  assistance  of  radar  and  each  aircraft  is  navigating  solely 
off  of  its  unaided  INS.  Assuming  the  two  aircraft’s  INS  errors 
are  independent,  Gaussian,  zero  mean,  and  have  a  1o  error  rate 
of  0.55  NM/hr,  y(0)  is  also  zero  mean,  Gaussian  with 


ayo  -  <°-55  taR)2  +  (0.55  taT)2 


Here  t  ^  and  t^  are  the  airborne  times,  in  hours,  of  the  receiver 


and  tanker.  For  an  airborne  time  of  10  hours  for  the  C-141  and 


5  hours  for  the  KC-135.  o^Q  is  computed  to  by  6.15  NM.  Of 


course  this  much  total  error  would  be  extreme.  Normally,  one 
or  both  aircraft  would  have  the  opportunity  to  update  its  INS 
somewhere  enroute,  or  the  rendezvous  track  may  be  defined  by 
radio  navigation  aids.  This  ’’worst  case"  is  included  to  get  a 
feel  for  the  limitations,  if  any,  of  the  proposed  system. 

The  case  in  which  radar  is  used  to  help  estimate  the 
receiver’s  position  is  also  modelled.  In  this  case 


E{y(0) }  =  yT(0)  +  AyQ 

where  y,p  is  the  tanker’s  y  coordinate  and  AyQ  is  a  manually 
inserted  estimate  of  the  lateral  distance  between  the  two 
aircraft.  The  accuracy  of  this  estimate  varies  with  range. 


Error  due  to  radar  azimuth  misalignment  decreases  linearly  with 
range.  Also,  at  shorter  ranges  the  radar  display  can  be  set  to 


a  smaller  scale.  No  data  on  the  accuracy  of  such  an  estimate 


could  be  found,  so  a  a  value  was  assumed.  From  personal 

V 

experience,  it  seemed  reasonable  to  say  that  a  crewmember  can 
estimate  bearing  from  a  radar  display  to  a  1  a  accuracy  of  2 
degrees.  At  60  NM  range,  a  1  degree  azimuth  error  results  in 
about  a  1  NM  cross-track  error.  The  relationship  between  azimuth 
and  cross- track  error  is  nearly  linear,  so  a  1  a  azimuth  error 
is  assumed  to  cause  a  a^Q  value  of  2  NM  for  a  60  NM  radar  range, 
and  a  a  value  of  1  NM  for  a  30  NM  radar  range. 

J 

Along-Track  Error 

The  receiver’s  along-track  position  and  velocity 
components  are  modelled  by 
x  =  v. 


■vD  cos  5  +  v  +  v 
R  wx  xe 


The  tanker  can  convert  the  receiver’s  planned  indicated 

airspeed  to  the  true  rendezvous  airspeed,  v^,  using  altitude 

and  air  temperature  data  from  its  central  air  data  computer  (CADC). 

INS  and  CADC  data  can  also  be  combined  to  compute  the  wind  velocity 

components,  v  and  v  .  Using  Eq(2),  5  is  then  calculated, 
wx  wy 

These  estimated  values  for  vD,  5,  and  v  are  only  computed  at 

ri  WX 

the  start  of  the  rendezvous  and  remain  constant.  In  Eq(4)»  v xg 
is  then  treated  as  the  only  random  variable.  v0 ,  5,  and  v 
could  be  continuously  updated,  but  then  the  spatial  correlation 
of  those  computed  values  and  v  would  have  to  be  accounted  for 


I 


in  the  truth  and  filter  models.  Certainly  this  would  be  an  area 
for  future  consideration,  but  it  was  considered  beyond  the  scope 
of  this  study. 

v  is  modelled  as  the  sum  of  a  random  constant  and  an 
xe 

exponentially- correlated  component.  The  constant  component  is 
primarily  due  to  instrumentation  and  pilot  bias  errors.  This 
bias  would  include  factors  such  as  tanker  errors  in  computing 
the  wind  components,  error  in  the  receiver’s  displayed  airspeed, 
the  receiver  pilot's  inability  to  maintain  the  correct  airspeed 
on  average,  and  differences  in  temperature  between  the  two 
aircrafts’  positions  when  v^  is  computed.  To  get  a  rough  idea 
of  the  magnitude  of  this  bias  error,  the  individual  error 
sources  were  assigned  variances.  The  figures  used  were  not 
explicitly  listed  as  standard  deviations;  approximate  la  valued 
had  to  be  surmised  from  published  CEP  values  and  tolerances, 
computer  flight  plans,  and  personal  experience.  The  following 
la  error  values  were  used:  tanker  ground  speed  -  1  knot  (Ref  1: 

6- 6),  tanker  and  receiver  true  airspeed  -  2  knots  each  (Ref  1: 

7- 5) »  airspeed  change  due  to  change  in  temperature  between  IP 

and  CP  -  1  knot  (Ref  3:  10-17,  30),  and  pilot  error  -  5  knots. 

Since  the  total  error  is  the  sum  of  each  of  these 

components,  and  they  are  each  assumed  to  be  independent,  the 

variance  of  the  total  error  is  the  sum  of  the  individual  variances. 

The  result  is  a  variance  of  35NM2  /hr2.  The  random  constant 

component  of  error,  v  ,  is  then  modelled  as  Gaussian,  zero  mean, 

xc 

with  a  standard  deviation  of  6  NM/hr. 
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The  correlated  component  is  due  mainly  to  wind  fluctua¬ 
tions.  Wind  data  is  difficult  to  interpret.  The  models  can  be 
very  complicated  and  vary  considerably  with  the  altitude  and 
frequency  ranges  considered.  A  basic  approximation,  though, 
is  that  wind  velocity  is  exponentially,  spatially  correlated 
(Ref  4:  2).  Aircraft  ground  speed  is  mainly  affected  by  just 

the  slowly  varying  winds  since  high  frequency  gusts  have 
no  net  effect  on  aircraft  displacement.  An  exponentially  corre¬ 
lated  wind  model,  with  a  correlation  distance  of  50  NM,  was 
adopted  from  (  Ref  4:  2). 

In  the  time  domain,  vx^»  the  fluctuating  component 
of  along-track  velocity  error,  is  modelled  by  the  equation 

Jxf(t)  =  -  —  vxf(t)  *  “f(t) 

where  T^  is  the  time  constant,  and  w^(t)  is  a  continuous  white 
noise  process.  T^.  is  found  from  the  equation 

rp  _  Correlation  distance 
Airspeed 

to  be  455  sec  for  a  C-141  traveling  at  396  knots.  Ala  value  of 

10  knots  was  chosen  for  vxf(t).  This  figure  seemed  reasonable 

after  examining  forecast  wind  variations  along  refueling  tracks 

on  actual  computer  generated  flight  plans  (Ref  3:  10-17,  30). 

The  appropriate  strength  of  the  zero  mean,  white,  Gaussian  noise 

w^.(t)  is  found  from  the  equation 
2a2 


q 


(Ref  5:  178) 


(5 


3 

I 

9 

t 

* 


4 


to  be  3.39  x  1 0”  ®NM2/ sec 3 .  The  resulting  vxf>(t)  process  is  zero 
mean,  as  it  should  be  if  the  wind  velocity  at  the  tanker’s 
position  is  considered  the  mean  wind.  Remember  that  v  ^  is  not 
the  total  wind,  it  is  just  the  fluctuating  component  of 
deviation  from  an  estimated  wind,  v 

wx 

Altitude  Separation 

In  the  truth  model,  position  and  velocity  components 
are  measured  in  the  local-level  coordinate  frame  discussed 
earlier.  However,  the  receiver  will  actually  be  displaced 
vertically  from  the  reference  plane  due  to  a  planned  altitude 
separation,  altimeter  errors,  and  the  curvature  of  the  earth. 

The  planned  altitude  separation,  aQ ,  is  usually  1000 
ft  or  0.164  NM.  Altimeter  error,  though,  will  cause  some 
deviation  from  aQ ,  even  if  both  aircraft  are  using  the  same 
altimeter  setting.  When  flying  at  a  constant  altitude  over  a 
relatively  short  distance,  the  error  is  essentially  constant 
(Ref  6:  4-4).  A  zero  mean,  Gaussian  random  bias  is  then  used 

to  model  the  net  error.  Maximum  acceptable  error  before  takeoff 
is  75  feet,  but  the  error  increases  slightly  with  altitude. 

For  the  simulated  rendezvous  altitude  of  25000  ft,  a  standard 
deviation  of  150  ft  or  0.0247  NM  is  assumed  (Ref  7:  1-5,  1-6). 

The  further  away  the  receiver  is  from  the  tanker,  the 
more  it  will  fall  below  the  tanker’s  local-horizon  reference 
plane.  The  amount  of  vertical  displacement  is  r  (1  -  cos  4>) 
where  rg  is  the  radius  of  the  earth,  and  <t>  is  the  angle  of 

20 
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central  arc  which  separates  the  two  aircraft.  The  angle  <J>  can  be 
found  from 

s 

d>  =  - 

r 

e 

where  s  is  the  arc  length  between  the  two  aircraft.  For  the 

short  distances  involved  -  less  than  100  NM,  s  is  approximated 

by  d,  the  straight  line  distance  between  the  tanker  and  receiver 

A  value  of  3444  NM  is  used  in  the  simulations  for  rg. 

The  altitude  separation  model  is  then 

a  =  a  +  a  +  a 
o  e  c 

where  a  is  the  total  altitude  separation,  aQ  is  the  planned 
separation,  a  is  altimeter  error,  and  a  is  vertical  displace- 
ment  due  to  the  curvature  of  the  earth. 

Incidentally,  the  velocity  component  v  ,  as  viewed  from 

A 

the  tanker,  will  also  be  affected  slightly  due  to  the  curvature 
of  the  earth.  The  receiver's  velocity  vector  is  oriented  out 
of  the  horizontal  plane  by  the  angle  <fc  described  above.  However 
the  maximum  horizontal  velocity  change  due  to  this  effect  is 
less  than  0.2  NM/hr  for  a  receiver  true  airspeed  of  400  NM/hr. 

At  the  initial  TACAN  ranges,  where  this  effect  is  the  greatest, 
the  final  estimator  design  only  achieved  ala  accuracy  of  about 
10  NM/hr,  so  ignoring  this  effect  in  the  truth  model  does  not 
affect  the  simulation  results. 


TACAN  Error 

With  air-to-air  TACAN,  the  tanker  can  measure  the 
straight-line  distance  to  another  TACAN  equipped  aircraft,  TACANs 
can  also  provide  a  bearing  measurement  to  TACAN  ground  stations, 
but  the  KC-135's  current  system  does  not  provide  bearing  capability 
for  air-to-air  operation  (Ref  10:  1).  The  TACAN  error  model 

used  in  (Ref  4:  3)  and  (Ref  8:  18)  was  adopted.  TACAN 

measurements  are  modelled  by  the  equations 

dT(t)  =  d(t)  +  dQ(t) 

and 

dg(t)  =  b  +  n(t) 

where  dT  is  the  TACAN  measured  distance,  d  the  true  distance, 
and  d  the  TACAN  error.  The  error  is  the  sum  of  a  random  bias, 
b,  and  a  short  correlation  time  component,  n. 

The  KC-135  was  recently  equipped  with  an  improved 
system,  the  AN/ARN  118  (V)  TACAN  made  by  Collins.  Test  results 
(Ref  9:  1),  published  tolerances  (Ref  10:  9),  and  flight 

test  results  of  a  similar  system  (Ref  11:  16)  indicate  the 

model  described  above  is  reasonable,  but  that  the  standard 
deviations  for  b  and  n  are  less  than  those  used  previously. 

The  bias,  b,  is  modelled  as  zero  mean,  normally  dis¬ 
tributed  (Ref  :  3),  and  having  a  standard  deviation  of  0.1  NM 

(Ref  10:  9).  n(t)  is  modelled  as  a  zero  mean,  Gauss-Markov 

process  (Ref  U'  3)  generated  by  the  system 

n(t)  =  -  n(t)  +  wn(t) 
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The  correlation  time,  T  ,  is  3*6  sec  (Ref  i:  3).  The  standard 
deviation  of  n,  for  the  new  TACAN,  is  0.016  NM  (Ref  9:  1) 

(Ref  11:  16).  The  appropriate  strength  of  the  continuous,  white 
Gaussian  noise  process  is  found  from  Eq(5)  to  be  1.422  x  10" “NM2/ 
sec . 

Summary 

The  entire  truth  model  is  summarized  below.  Random 
variables  are  underscored  with  a  tilde,  and  mean  and  standard 
deviation  values  are  denoted  by  m  and  a  respectively.  Non- 
random  variables,  such  as  vR,  are  determined  from  planned  rendez¬ 
vous  parameters,  and  INS  and  CADC  data.  The  listed  values  are 
those  used  in  the  simulations. 


6  = 


sin 


6  =  Drift  correction 

vwy  =  y  component  of  wind 

vR  =  Receiver  true  airspeed 

vR  =  0.110  NM/sec  (396  NM/hr) 

£(t)  =  Yyc  +  ?y(t) 

^(t)  =  Receiver's  y  coordinate 


(6) 


(7) 


No  radar  aiding: 


With  radar  aiding: 


myo  =  M0)  +  Ay0 

where  =  Tanker  y  coordinate 

AyQ  =  Estimated  offset 

Oy0  =  1  NM  for  30  NM  radar  range 

Constant  component  of  lateral  drift 

mvyC  =  0,  avyc  =  2.161  x  10"“NM/sec  (0.778  NM/hr) 

Continuous  white  noise 

m  =  0,  q  =1.13x10"  5NKV sec 
wy  My 

x(t)  =  -vR  cos  6  +  vux  +  vxc  *  vxf(t) 

Receiver  x  coordinate 
x  component  of  wind 


Constant  component  of  x  velocity  deviation 


m  =  0 ,  a  =  2.78  x  10“3NM/sec  (10  NM/hr) 
vxc  uvxc  ' 


Fluctuating  component  of  deviation 


m vxf  =  °'°vxf  =  1,67  x  10~3NM/sec  (6  NM/hr) 


xf 


(t)  = 


vxfU)  +  wf(t) 


Correlation  time 
T -  =  455  sec 


Continuous  white  noise 


wf 


=  0,  qf  =  3.39  x  1 0~  ®NM2/ sec  3 


a  =  +  a„ 

~  o  ~  e  c 

Altitude  separation 
Planned  separation 

aQ  =  0.164  NM  (1000  ft) 

Altitude  error 

=  °»  °  =  0.0274  NM  (150  ft) 

Separation  due  to  curvature  of  the  earth 

a  =r  (1-cos(d/r)} 
c  6  6 

3444  NM  ' 

Distance  between  aircraft 
(t)  =  b  +  n(t) 

5  ^ 

TACAN  measurement  error 
TACAN  bias 
rab  =  0*  ab  =  °.1 

Short  correlation  component  of  error 

m  =  0,  a  =  0.0160  NM 
n  n 

t)  =  -  — — n(t)  +  w  (t)  (12) 


i  a* 


As  stated  earlier,  the  linear  form  of  this  resulting 
dynamic  equation  motivated  the  specific  choice  of  coordinate 
system . 

The  discrete  measurements  at  time  t^  are  described  by 
T  ( )  =  J{x(  tT)  -  xT(t.)}2+  (^(ti)-yt(ti)  }2+  {aQ  +  ag  +  aQ}  2 

+  b  +  n(t.)  (U) 

With  the  truth  model  for  state  dynamics  and  measurements 
defined,  attention  will  now  be  turned  to  designing  an  estimator 
to  give  the  necessary  navigational  inputs  to  the  rendezvous 
control  algorithm. 
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IV  Estimation 


The  design  goal  is  a  system  which  would  enable  the  tanker 
to  rendezvous  more  accurately  with  the  receiver.  The  better  the 
tanker  estimates  the  receiver’s  position  and  velocity,  the  better 
he  can  determine  where  and  when  to  turn.  Also,  a  continuously 
updated  estimate  of  the  receiver’ s  state  would  enable  the  tanker 
to  use  closed  loop  control  during  the  final  turn  -  a  capability 
which  does  not  exist  with  the  present  method.  To  provide  such 
an  estimate,  a  reduced  order,  extended  Kalman  filter  was  designed 
Extended  Kalman  Filter 

The  descrete,  linear  Kalman  filter  (Ref  5:  Ch  5)  is  a 
commonly  used  estimation  method.  Its  estimate  of  a  state  vector 
x  is  based  on  an  assumed,  linear,  stochastic  dynamics  model  of 
the  form 

x(t)  =  F(t)x(t)  +  B(t)u(t)  +  G(t)w(t)  (15 

where  u(t)isa  deterministic  forcing  function  and  w(t)  is  a 
continuous  white  noise,  and  a  linear  measurement  equation 

l(t1)  =  H^)  x( t± )  +  v^)  (16 

where  _z(  t^ )  is  a  discrete  measurement  vector,  and  v(t^)  is  a 
discrete  white  noise.  The  Kalman  filter  also  assumes  x(tQ), 
w(t),  and  v(t^)  to  be  Gaussian  and  mutually  independent. 

The  dynamic  and  measurement  equations  developed  for  the 
truth  model,  Eqs(13)  and  (14)>  Tit  the  standard  Kalman  filter 
form  and  assumptions  except  the  measurement  equation  is  nonlinear 


The  extended  Kalman  filter  is  one  means  of  handling  such  non- 
linearities  (Ref  12:  Sec  9.5).  It  was  chosen  over  higher  order 
nonlinear  filters  (Ref  12:  Ch  12)  for  this  initial  study  because 
it  is  more  likely  to  be  operationally  employed  since  it  poses  less 
computational  burden. 

In  the  extended  Kalman  filter,  the  nonlinear  dynamic  and 
measurement  equations  are  linearized  about  the  most  recent  state 
estimate.  The  result  is  a  linear,  state  and  measurement  pertur¬ 
bation  model.  Now  the  standard  Kalman  filter  update  equations 
can  be  used  to  estimate  the  perturbation  state  5  x>  and  the  state 
vector  x  can  be  updated  by 

x(tt)  =  x(t7)  +  5x(tt) 

where  the  hat  notation  indicates  an  estimate,  and  the  time 
arguments  tT  and  tt  distinguish  between  state  estimates  before 
and  after  update. 

For  a  problem  such  as  this  one,  in  which  only  the 
measurement  equation  is  nonlinear,  the  extended  Kalman  filter 
can  be  summarized  as  follows: 

Dynamics  Model: 

$(t)  =  F(t)x(t)  +  B(t)u(t)  +  G(t)w(t)  (17) 


Measurement  Model: 


Statistical  Description  of  Uncertainties: 

E{x(to)}  =  xQ  E{(x(to)-x0)(x(to)-xo)  }  =  PQ 

E{w(t)}  =  0  E{w(t)w^  (t+x)}  =  ^(t)  6(x) 

E{v(ti)}  =  0  E{v(ti)vT(ti)}  =  R(ti) 

E{ v(ti)wT(t)}=  0 

where  5^  is  the  Kronecker  delta  function. 

State  Propagation  Between  Updates: 
x(tT)  =  J(ti,ti_1)  x(tt_x) 


+J  1  £  (titx)B(x)u(x)dx 


(19) 


p(t")  =  iCti.t^jjpet^j)  gT(ti,ti_l) 


+ /  i(ti,x)G(x)<l(x)GT(x)£T(ti,x)dx 


Vi 


where^is  the  state  transition  matrix  defined  by 

ft4(t,t0).  £(t)  K‘’V'  £<W  *  - 

Measurement  Update: 

K(t.)  ={P(t‘)HT(t. )  H(t.)P(t‘)HT(t.)  +  R(t.)} 


where 


H(t.)  = 


3  h  (x,  t..) 


x  =  x(tT ) 


(22) 
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(23) 


Sx(tt)  =  K(ti){  z{  ti)-h  (  x(tT ) ,  tJ  } 


x(tt) 

=  x(t7)  +  5x(tt) 

(24) 

p(t*) 

=  P(t')  -  K(ti)H(ti)P(tr) 

(25) 

Dynamics  Model 

The  truth  model  equations,  (13)  and  (14),  could  be  used 
for  the  filter  model  as  well.  Then  filter  design  would  simply 
be  a  straightforward  application  of  equations  (17)  through  (25). 
However,  unnecessary  states  mean  unnecessary  computations,  and 
a  trade-off  should  be  made  of  computer  loading  versus  filter 
performance.  An  analysis  of  the  error  states  in  the  truth  model 
indicates  that  the  filter  can  probably  be  reduced  to  just  three 
states  without  a  significant  loss  of  accuracy. 

First  of  all,  the  altitude  model  can  be  simplified. 

The  amount  of  altitude  separation  due  to  the  curvature  of  the 
earth  can  be  large,  but  it  has  little  effect  when  estimating 
velocity  or  position  in  a  horizontal  plane.  For  example,  for 
a  distance  between  aircraft  of  100  NM,  a  equals  1,45  NM. 

However,  ignoring  a  would  only  c*">.se  a  velocity  estimate  error 
of  0.17  knots  and  a  position  estimate  error  of  65  ft;  at  a  closer 
range  of  30  NM,  the  errors  would  only  be  0.01  knots  and  1.73  ft. 

A  similar  analysis  of  the  effect  of  altimeter  error  shows  that 
it,  too,  can  be  ignored.  Even  at  the  nominal  rendezvous  completion 
range  of  2  NM,  a  3o  altimeter  error  of  450  ft  would  only  change 
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the  horizontal  position  estimate  by  45  ft.  So,  only  the  planned 


altitude  separation  term,  a  ,  is  retained  in  the  filter  model. 

The  filter  model  was  further  simplified  by  dropping  the 
measurement  noise  components  from  the  state  vector.  The  correlated 
component,  n,  with  a  la  value  of  0.016  NM,  is  small  compared  to 
the  bias,  b,  which  has  ala  value  of  0.1  NM.  Also,  its  short 
correlation  time  makes  it  appear  to  the  filter  as  a  discrete, 
uncorrelated  noise.  For  example,  with  a  sampling  time  interval 
of  10  seconds,  the  covariance  of  2  consecutive  samples  of  n  is 
only  1.59  x  10"5NM2.  It  would  seem  advantageous  to  be  able  to 
estimate  the  bias.  This  was  attempted;  but  as  will  be  discussed 
further  in  Section  VI,  simulations  indicated  that  b  is  relatively 
unobservable.  When  it  is  retained  as  part  of  the  state  vector, 
its  variance  decreases  very  little  and  no  noticeable  improvement 
in  filter  performance  is  attained.  Therefore,  in  the  final  filter 
design,  no  attempt  is  made  to  estimate  explicitly  b  or  n. 

In  the  truth  model,  the  along-track  states  are  modelled 


x(t)  = 


-v^  cos  6  + 


v  + 
xc 


vxf(t) 


But  it  would  be  difficult  for  an  estimator  to  differentiate 


between  the  random  constant  vxC  and  the  slowly  varying  One 

technique  is  to  replace  a  bias  plus  long  correlation  time  error 


with  a  random  walk  model  (Ref  1  3 1  2-3). 


vxc  +  vxf  therefore 


replaced  by  a  single  state,  v  ,  which  is  modelled  by 

X  G 


where  wx  is  a  continuous  white  noise. 

Finally,  the  cross-track  model  is  simplified.  The  truth 
model  equation  is 


y(t)  =  vyc  +  wy(t)  (7) 

However,  the  magnitude  of  the  random  constant  vyc  is  so  small 
that  no  attempt  is  made  to  estimate  it  explicitly.  Instead,  it 
is  dropped  from  the  filter  state  vector  and  the  strength  of  the 
white  noise  is  increased. 


What  is  left  is  a  simple,  three  state  dynamics  equation 


which  only  models  the  dominant  states: 


0  0 

1  0 

0  1 

G  w 

These  states  are  also  the  only  input  data  required  by  the  controller. 

In  the  reduced  order  filter  model,  there  is  no  exact  way 
of  setting  the  noise  strengths.  Experimentation,  or  tuning,  will 
determine  which  values  are  chosen.  Some  reasoning  can  be  used, 
though,  to  arrive  at  initial  trial  values  for  £. 
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When  replacing  one  model  with  another  simpler  one,  it 
would  seem  appropriate  for  the  outputs  of  both  models  to  be 
changing  value  in  some  similar  manner.  For  an  exponentially 
correlated  component  replaced  by  a  random  walk,  an  empirically 
found,  useful  technique  is  to  set 


where  a2  is  the  correlated  process  variance  and  t  is  the  correla 
tion  time  (Ref  13:  2-3).  This  method  was  used  to  find  an 

appropriate  strength  for  wx: 

o2  r. 

_  - vxi —  _  1#69  x  i  o~  8NM2/ sec  3 


A  similar  insight  was  used  to  set  q  .  The  one  sigma 

«/ 

cross- track  position  error  was  said  earlier  to  be  0.55t  NM  wher 

a 

t  is  the  time  in  hours  since  alignment.  The  net  cross-track 

cL 

error  at  rendezvous,  then,  has  a  variance  of 

a2y  (t)  =  (0.55taT)2  +  (0.55taR)2=  0.3025(taT2+  t&R2) 


where  ta^  and  t&R  are  the  airborne  times  of  the  tanker  and 
receiver  at  rendezvous.  By  taking  the  derivative  of  a2  (t)  with 

y 

respect  to  time,  it  is  found  that  the  variance  at  rendezvous  is 
changing  at  the  rate  of  0.605  ^aT  +  W  NM  /hr.  Because  the 
variance  of  the  random  walk  filter  model  is  q  (t-tQ),  a 
value  of  qy  also  equal  to  0.605  ( ta^»  +  t&R)  would  result  in  a 
filter  model  changing  variance  at  about  the  same  rate. 
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For  a  relatively  long,  combined  airborne  time  at  rendezvous  of 
15  hours,  q  =  2.52  x  10-3NM2/sec. 

y 

Measurement  Model 

For  the  tree  state  filter  model,  Eq(14)  simplifies  to 
dT(ti)  =  v  {  x(  tj_ )  -  xT(ti)}2  +  {  ytt^.  )-yT(ti)}  2+  aQ2+v(ti)  (27) 
z(ti)  =  h(x(ti),  t^  +v(ti) 

Since  this  is  a  simplification  of  the  truth  model,  a 
suitable  variance  for  the  white  noise,  v,  will  also  have  to  be 
found  through  trial  and  error.  Primarily,  v  is  replacing  b  +  n, 
so  a  reasonable  first  guess  is  to  set 

R  =  °  b + 

For  the  variances  of  b  and  n  used  in  the  truth  model,  the  initial 
filter  value  of  R  is  calculated  to  be  0.101  NM2 

Initiali zation 

The  tanker  starts  the  filter  initialization  process  by 
estimating  vR,  vyx,  and  v ^  -  the  receiver’s  true  airspeed  and 
the  components  of  wind  affecting  its  ground  velocity.  These 
would  be  calculated  from  the  planned  indicated  airspeed  for  the 
receiver;  INS  computed  ground  speed,  track,  and  heading;  and 
CADC  measured  altitude  and  air  temperature.  The  receiver's 
drift  is  then  estimated  from  Eq(2).  Or,  if  the  receiver  has 
readouts  of  any  of  those  values  onboard,  they  can  be  relayed 


over  the  radio  and  inserted  manually  into  the  computer.  The 

best  guess  of  the  receiver’s  x  velocity,  then,  is  -vDcos  5 +  v  ; 

rt  wx 

and  estimated  deviation  from  this  value,  v  (0),  is  zero. 

The  tanker  then  calculates  and  attempts  to  fly  the  correct 

offset.  If  no  radar  aiding  is  used,  y(0)  =  0.  If  an  estimate 

of  the  lateral  distance  between  the  two  aircraft,  AyQ  ,  is  made 

with  the  help  of  radar  or  other  means,  then  it  is  manually 

inserted  and  y(°)  =  yT(0)+  Ay0» 

The  initial  covariance  matrix  values  are  also  tuned  for 

best  performance.  Since  P  is  symmetric,  only  upper  triangular 

elements  are  calculated.  The  initial  variance  for  y  can  be 

taken  directly  from  the  truth  model:  P33(0)  =  a2  .  Also, 

because  v  =  v  +  v  it  would  make  sense  to  choose  P22(0)  = 
xe  xc  xf 

°vxc  +  °vxf*  The  tanker  would  initially  have  very  little  infor¬ 
mation  about  the  correlation  of  x,y,  and  v  ,  so  Pi2(0)  and  P23(0) 

X  G 

are  set  to  zero. 

The  only  initial  conditions  not  yet  specified  are 
x(0),  Pn(0),  and  Pi3(0).  A  TACAN  measurement  at  t=0  can  be 
used  to  approximate  those  values.  Approximate  mean  and  covariance 
relations  are  derived  using  a  Taylor  series  expansion  about 
y(0)  and  the  initial  TACAN  measurement. 

The  actual  distance,  d,  is  related  to  x,  y,  and  a  by 

d  =  J  ( x-  Xip )  +  (y-yT)z+  a2 


The  effects  of  altitude  deviations  from  aQ  are  ignored,  as 
discussed  earlier,  and  the  equation  above  is  rearranged  to  form 
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a  function  Ax(d,Ay): 

Ax(d,Ay)  =  J  d2  -  Ay*  -  aj"  (2? 

where  Ax  =  x-xT  and  Ay  =  y-yT>  Now  Eq(27)  is  expanded  to  first 
order  about  AxQ  =  Ax( d^, , Ay)  where  d^  is  the  TACAN  measured 
distance,  and  Ay  =  y-y^: 


Ax(d,Ay)  =  Ax  +  — ^ —  (d-dT)  + 

°  3d  1 


—  (Ay-Ay) 


d  =  da 


Ay  =  Ay 


d  =  d„ 


Ay  =  Ay 


d-drp  is  the  TACAN  error,  dg;  and  Ay  -  Ay  =  ( y- yT ) -  ( y-yT )  =  y-y 
which  is  the  cross-track  error  y  .  The  derivatives  are  evaluated 
and  Eq(29)  reduces  to 

Ax( d, Ay )  =  ix0  +  dT  d  -  Ay  y 


Now  x(0),  Pix(O),  and  Pi  3(0)  can  be  approximated  to  first 
order  by  taking  the  appropriate  expected  values  and  making  use 
of  Eq(30).  Remember  that  d,p  and  Axq  are  not  random,  and  that 
dg  and  yg  are  assumed  to  be  zero  mean  and  independent. 


x(0) 


X  =  Xrp  +  Ax 


=  X,p  +  AXq  +  drj 


de  -ii. 
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1 


r 


0 


f 


1  (t..T) 


0  0 


0  qx  0  (t.-T)  1 


10  0  q 


Since  P  is  symmetric,  only  the  upper  triangular  elements  need  to 
be  calculated.  The  matrix  equation  above  reduces  to  the 
following  relations: 

Pii  =  P  i,  i  +  2  Pi2^t^  +  P22  At?  +  -y  q^  At?  (3 

P12  =  P12  t  P 2 2 A t .  +  £  q  At?  (3 

P22  =  P22  +  qx  At^  (3 

P 1 3  =  P13+P23  At^  (3 


=  P  2  : 


=  P 


33  +  q  At^ 


where  the  P  terms  on  the  right  are  the  upper  triangular  elements 
of  P(tt_1)>  and  the  P~  terms  are  elements  of  P(t7). 


Update 


From  the  filter  measurement  model,  Eq(27), 


h^x^)}  =  {[xCt^  -  xT(t.  )]  2  +  [y ( t± )  -  yT(ti)] 


Ax(  t .  )  2  +  Ay(ti)2  +  a 


Then,  from  Eq(26), 


3h  {  x(  t^  )} 

3  x 

A 

X  =  X 

Ax(ti)  q  Ay(ti) 

W 

df(ti 

(38) 


where  is  the  filter  computed  distance 

df(ti)  =  Jax(  t^ )  2  +  Ay^)2  +  a”2 

The  state  and  covariance  updates  are  now  just  a  straight¬ 
forward  application  of  equations  (21)  through  (25).  The  result¬ 
ing  update  relations  are  listed  in  the  next  section. 

Sometimes  when  measurements  are  nonlinear,  filter  perfor¬ 
mance  can  be  improved  by  adding  a  bias  correction  term  to  the 
predicted  measurement  value,  h(x)  (Ref  12:  225).  For  a  single 
measurement,  the  appropriate  bias  correction  is  given  by 

r  ... 


bc  =  *  tr 


32h(x) 

3x2 


The  correction  term  for  this  problem  is  then  found  to  be 

1  Ax2  1  Ay2 

b  =  4(0 — - )Pn+  ( —  -  -  )  P33} 

c  d„  d  3  d „ 


df  df3 


Now  the  predicted  measurement  value  for  forming  the  residual  in 
Eq(22)  is 

z(ti)  =  h(x(ti)}  +  bc('ti) 
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In  effect,  z  was  just  estimated  by  expanding  h(x)to  second  order 
and  taking  the  expected  value.  Eq(22)  now  becomes 

5x  (tt)=  K(  (z(ti)  -  z(ti)  }  (41 

Estimation  Algorithm 

The  entire  estimation  algorithm  is  now  summarized. 
Anytime  prior  to  rendezvous,  CP  coordinates,  rendezvous  true 
course,  receiver  indicated  airspeed,  and  aQ  are  loaded  into 
tanker's  navigation  computer.  At  each  measurement  time,  INS 
computed  tanker  position  is  transformed  to  the  coordinate  system 
illustrated  in  Figure  2  on  page  10. 

At  start  of  rendezvous: 

1.  INS  and  CADC  data  are  used  to  compute  vD ,  v  ,  v  ,  and 

r  wx  wy 

6  =  sin 


2.  u  =  -vR  cos  6  t  vux 


3.  Tanker  computes  and  attempts  to  establish  correct  offset. 

Initialization: 

1.  If  an  estimate  of  the  receiver's  actual  offset  can  be 
obtained  from  radar  or  other  means,  Ay(0)  is  inserted 
manually.  If  no  update  is  made,  Ay(0)  equals  desired  offset. 

2.  The  filter  is  started  with  an  initial  TACAN  measurement, 

dT(0). 

3.  x  and  P  are  initialized  as  follows.  The  time  argument  for 


all  variables  is  t=0. 


Ax  =  /  dT2  -Ay2  -  a 


xT  +Ax 


yT  +  Ay 


— —  d£  (a,  2  +  a  2)  +  Ay2  a  2 
^2  [_  T  b  n  *  yo 


a  2  +  a  2 
vxc  vxf 


Av  a  2 
— 2 —  yo 

Ax 


Propagation: 

x  and  P  are  propagated  from  sample  time  t^  to  t^  with  the 

following  equations.  Minus  superscripts  indicate  tT  values. 

the  variables  on  the  right  side  of  the  equations  are  tt  values, 

At  =  t.  -  t. 

l  l- 1 

x"  =  x  +  (v  +  u)At. 

Xc  X 


V  =  V 

xe  xe 


y  =  y 
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Pii"  =  Pn  +  2Pi2Ati  +  P22At|  +  -  ^xAti3 

Pl2~  =  P i 2  +  P  2  2  A t .  +  i  q  At.2 
P  2  2  ~  =  P22  +  ^xAti 

P 1  3  =  P 1  3  t  P  2 3At^ 

P  2  3  =  P  2  3 

P33  =  P 3  3  t  q  At^ 

Update: 

At  each  sample  time,  a  TACAN  measurement  d^,  is  used  to  update 
x  and  P.  The  following  relations  are  derived  from  equations 
(21)  through  (25)  with  a  bias  correction  term  added.  x  and  P 
elements  on  the  right  side  of  the  equations  are  t7  values. 

®  Ax  =  x  -  xT 

A  A  A 

Ay  =  y  -  yT 

d^  ='1ax2  +  Ay2  +  aQ2 

Hi  =  Ax/d^ 

H2  =0 

H3  =  Ay/df 

A  =  H 1 2P i 1  +  2H1H3P13  +  H32P33+  R 


PH  1 

-  H iP  1 1 

+ 

H  3P 1 3 

PH  2 

=  HiP12 

+ 

H3  P23 

PH3 

=  H iP  1  3 

+ 

H  3P  3  3 

Ki 

=  PHi/A 
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*<> 


9) 


k2 

k3 


z 


ph2/a 

ph3/a 

i  {(1/df  -  Ax2/df3)Pi!  +  d/df  -  Af/df3)P33} 

df  +  bc 

dT  -  z 

x  +  Kir 

v  +  K2r 
xe 

y  +  K  3r 
Pn-KxPHi 
Pi2-KiPH2 
P  2  “K  2  PH  2 

P i 3-K iPH3 

p23-k2ph3 

p33-k3ph3 

This  algorithm  is  used  to  continually  propagate  and 


update  the  receiver' s  state  estimate  for  as  long  as  needed  by 
the  controller. 


V  Control 


Algorithms  were  developed  to  calculate  the  correct  offset 
distance,  the  appropriate  time  for  the  tanker  to  start  its  turn, 
and  bank  angle  commands  during  the  turn.  For  this  initial 
feasibility  study,  the  control  laws  are  based  on  a  deterministic 
dynamics  model  with  the  current  receiver  state  estimate  as  input. 
Tanker  Dynamics 

Eq(1)  is  appropriate  for  deriving  the  kinematic  relations 

fox*  the  tanker  as  it  was  for  the  receiver: 

v  =  v  /  +  v 

—  —a/  w  — w 


where  v  is  velocity  relative  to  the  ground,  v  ,  is  the  velocity 
of  the  aircraft  with  respect  to  the  air  mass,  and  v^  is  the 
velocity  of  the  wind.  The  tanker's  heading-,  T  ,  is  measured 
counter  -  clockwise  from  the  positive  x  axis  as  shown  back  in 
Figure  2.  Then  for  a  tanker  true  airspeed  of  vT,  the  components 

of  -v„ are 

va/wx  ’  VT003 
va/wy  '  VT  sin'1' 


The  components  of  the  tanker’ s  velocity  in  the  pseudo 
fixed  frame  of  the  INS  are  therefore 

vTx  =  VT00S'1'  +  vwTx 

„Ty  =  vTsin  T  t  vwTy 

where  vwTx  and  vwT^  are  the  components  of  wind  at  the 
current  position. 


ground 


(4D 

(42) 

tanker’ s 
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A  simplified  free  body  diagram  of  a  KC-135  in  a  left  turn, 
with  positive  left  bank  angle  <t>,  is  shown  in  Figure  3.  L  is  lift, 
m  the  tanker's  mass,  and  g  the  acceleration  due  to  gravity.  In 
a  level  turn,  the  vertical  component  of  lift  balances  the  weight. 


mg,  so 


L  cos  <}>  -  mg  =  0 


(43) 


In  a  turn,  the  tanker  has  a  centripetal  acceleration,  ¥  v^,, 

which  is  caused  by  the  horizontal  component  of  lift.  The  resulting 


equation  of  motion  is 

L  sin  T  =  mTvrp  ( 44) 

Equations  (43)  and  (44)  are  now  combined  to  form  a  relation  for  $ : 
&  _  g  tan* 


Some  simplifications  were  employed  in  developing  this 
model,  but  these  equations  represent  the  principle  dynamic  effects 
throughout  the  turn.  The  air  mass  and  INS  coordinate  frame  are 
assumed  to  be  non-accelerating,  which  they  virtually  are  when 
compared  to  the  accelerations  of  a  turning  aircraft.  Also,  the 
velocity  vector  is  assumed  to  be  along  the  tanker's  longitudinal 
body  axis.  For  a  KC-135  at  cruise  airspeed  and  shallow  bank 
angles,  the  angle  of  attack  is  small;  and  its  effect  on  bank  and 
heading  calculations  is  assumed  to  be  negligible  for  the  purposes 
of  this  problem.  As  in  the  filter  design,  an  implementable 
control  law  has  to  be  based  on  a  model  of  just  the  dominant 


effects. 
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Turn  Point 


m 


The  geometry  of  the  tanker’s  turn  is  illustrated  in 
Figure  4»  For  a  constant  bank  angle,  constant  airspeed  turn,  the 
tanker's  flight  path  relative  to  the  air  mass  is  circular  with 


radius  of  turn 


( 46) 


The  tanker  wants  to  roll  out  on  the  refueling  track  with  zero 
lateral  drift.  Then,  as  the  receiver  did,  it  needs  to  apply  a 
drift  correction;  and  its  final  heading  will  be 


=  it  +  6 

where 


(47) 


From  any  initial  heading,  ¥  ,  the  change  in  the  tanker's  position 

during  the  turn  will  equal  the  change  in  position  within  the  air 
mass  plus  the  change  due  to  movement  of  the  air  mass  over  the 
ground.  Then 

AXiji  —  Xrj,  —  ^'j'o 

*  rT  [sin  (»  +  «)  -  8in»0]  +  vuIxT 

=  rT  [-sln  5-  sln  *„]  +  vwTxT  <48) 

where  and  x^Q  are  the  tanker  final  and  initial  x  coordinates 

and  T  is  the  duration  time  of  the  turn.  T  can  be  found  from 


T 


A  y 
♦ 
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7T  +  5-f  0 
Y 


(49) 


By  substituting  Eq(45)  into  Eq(46),  it  is  found  that 


g  tan  <t> 

Also,  substituting  Eq(45)  into  Eq(49)  results  in 

Vm 

T  =  - ±- 

g  tan  <J> 

Now  these  relations  for  r^,  and  T  are  used  in  Eq(48)  to  derive 


Axt  = 


g  tan<t> 


{  v^(- sin  6  -  sin  V Q) 


+  ("+5-  V  vwTx  ) 


(52) 


Similarly,  it  can  be  shown  that 


AyT  =  ~  yTc 


=  rT  [cos  6  +  cos  ¥0]  +  vwTyT 


g  tan  <J> 


{vT (cos  6  +  cos4'q)  +  (tt+5-  fQ)vwTy}  (54) 


Eq(54)  can  be  used  directly  to  determine  the  nominal 
offset  distance.  For  the  tanker  to  fly  towards  the  receiver 
parallel  to  the  refueling  track,  a  drift  correction  of  -6  must 
be  applied.  The  change  in  sign  is  because  6  is  calculated  for 
the  rollout  heading,  and  the  tanker  is  initially  flying  in  the 
opposite  direction.  Then  f  =  -6  and  the  desired  offset  is 


g  tan  <b 


{  2vtcos  6  +  (tt+  25)  v™  } 
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where  <t>  is  the  nominal  bank  angle. 

The  present  procedure  uses  a  nominal  bank  angle  of  30 
degrees.  But  30  degrees  is  also  the  tanker's  maximum  bank  angle 
under  normal  conditions.  Therefore,  a  shallower  nominal  bank 
angle  must  be  used  if  bank  angle  corrections  are  to  be  applied 
during  the  turn.  With  no  cross-wind,  a  nominal  bank  angle  of 
25  degrees  requires  about  a  9.6  NM  offset.  Nominal  offsets 
greater  than  that  were  considered  unreasonable  due  to  airspace 
restrictions  imposed  by  air  traffic  control. 

The  tanker  attempts  to  establish  the  nominal  offset 
distance,  but  the  computed  time  to  turn  should  be  based  on  the 
tanker's  actual  offset  and  heading  -  not  the  nominal.  For 
example,  the  tanker  may  actually  be  wider  than  the  computed 
offset  and  its  initial  heading  other  than  -6  while  it  is  correct¬ 
ing  to  course.  In  that  case,  it  should  turn  earlier  and  with 
a  shallower  bank  angle  than  nominal.  The  ability  to  compute 
other  than  nominal  solutions  does  not  exist  with  the  present 
method  but  could  be  implemented  with  the  proposed  system. 

To  compute  4>r,  the  required  bank  angle  to  roll  out  along 
the  receiver's  estimated  flight  track,  y  -  y^,  is  substituted 
for  Ayrp  and  Eq(54)  is  solved  for  <t>: 

$r  =  taRl  iTyryT){  VT(C°S  5  +  C0S  'l’o)  +  (lT  +  5"fo)  vwTy}  (56) 
Given  <t>r,  the  time  required  for  the  turn  is 
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T  = 
r 


(tt  +  6-4'o)  vt 
g  tanb 


which  is  derived  from  Eq(5l).  If  the  tanker  turns  at  the  current 
time,  its  x  coordinate  at  rollout  will  be 

xTf  =  XT  +  AxT 

where  Ax,p  is  found  from  Eq(52).  The  receiver's  position  at 
tanker  rollout  is  estimated  by 


A  A  .  A  rp 

xf  =  x  +  vx  Tr 


where 


vx  =  (-vRM3«t  VMX)  t  vxe  (58) 

The  tanker  wants  to  roll  out  2  NM  in  front  of  the  receiver, 
If  it  turns  too  early,  its  miss  distance  in  front  of  the  receiver, 

v  is 

xm  =  xf  -2  -xTf 

But  x  is  decreasing  at  the  rate  of  the  two  aircraft's  closing 
m 

speed,  vT  -v  .  The  estimated  time  to  go  until  the  turn  should 

XX  X 

be  started,  TTG,  is  then  given  by 


TTG  = 


vTx  "  Vx 


TTG  is  counted  down  between  updates,  and  the  tanker  turns  with 


bank  angle  br  when  TTG  £0. 
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The  algorithm  for  determining  the  turn  point  is  now  summarized: 
A.  The  tanker  computes  and  attempts  to  establish  the  correct 


offset  for  a  nominal  bank  angle  <J> 


1  x  -  Q4  =  i  wT" 

I .  o  =  sm  —  ■  '■ 


2.  os  = 


g  tan  4> 


{2v^*cos6  +  (tt+  26)vwrp^.} 


The  required  bank  angle  and  time  to  go  until  turn  are 
updated  based  on  actual  tanker  position  and  heading: 


1.  A4>  =  TT  +  S-'i* 


( 60a) 


2.  <J>  =  tan 

r 


g(y-yT) 


{  Vrp  (  cos5+cos  4^)  +  Al'v^  } 


(60b) 


3.  T 


AH'*  v„ 


g  tan  <t> 


( 60c) 


4-  •  A  Xrp  * 


g  tan 


{vT(-sin5-sin4'o)  +AT  •  vwTx> 


(60d) 


5  •  Xrp  £•  ~  Xfp  t  A  Xpp 


( 60e) 


6.  v  =  (-vdcos5+  v  )  +  v 


(60f) 


n  ^ 

7.  x«  = 


x  +  v  T 
x  r 


( 60g) 


8.  x_  = 


xf  -  2 


( 60h) 


9.TTG  = 


vTx  "  vx 


(60i) 
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Closed  Loop  Control 
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As  was  stated  previously,  one  of  the  primary  potential 
benefits  of  the  proposed  system  is  that  it  would  enable  the 
tanker  to  use  closed  loop  control  during  the  final  turn  to 
rendezvous . 

The  present  procedure  is  entirely  open  loop:  the  tanker 
maintains  a  constant  30  degrees  of  bank  until  reaching  the 
refueling  heading.  A  simple  improvement  would  be  to  apply 
cross-track  control  only.  At  each  update  time,  the  tanker  would 
compute  the  required  bank  angle  to  roll  out  in  front  of  the 
receiver  using  Eq(56).  The  result  should  be  decreased  cross¬ 
track  error  at  rollout,  which  would  mean  less  lateral  maneuvering 
would  be  required  of  the  receiver  to  complete  the  join-up. 

Using  bank  angle  control  to  minimize  the  along-track 
error  also,  is  more  complicated.  That  is  because,  in  general, 
the  tanker  cannot  get  from  its  present  position  and  heading  to 
a  zero  error  rollout  condition  with  a  constant  bank  angle.  A 
method  of  solving  the  problem  which  was  used  in  the  simulations, 
is  based  on  the  following  basic  algorithm: 

1.  Compute  <t>r  to  roll  out  on  the  receiver’s  estimated  track. 

2.  Determine  the  x  miss  distance  if  <j>  is  used. 

r 

3.  Adjust  4>r  based  on  x  miss.  If  the  tanker  would  roll  out 

in  front  of  the  receiver,  decrease  the  bank  angle  and  vice- 
versa. 

The  basic  idea,  then,  is  to  see  if  an  exact  solution  exists  from 
the  present  position  and  heading.  If  one  does  not,  adjust  the 
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bank  angle  to  fly  towards  one. 


The  first  two  steps  of  this  algorithm  are  solved  by 

equations  (60a)  through  (60h).  The  remaining  question  is  then 

how  much  to  correct  <J)  .  An  ad  hoc  method  was  used  based  on  the 

r 

following  insights.  If  the  tanker  would  roll  out  in  front  of 
the  receiver,  it  needs  to  fly  a  longer  path  to  the  rendezvous. 
Similarly,  if  it  would  roll  out  behind  the  receiver,  more  cuttoff 
and  a  more  direct  route  should  be  flown.  A  measure  of  the  effect 
of  a  bank  angle  change  at  the  present  position,  is  the  change 
in  arc  length  it  would  cause  if  held  until  rollout.  The  control 
law  was  then  based  on  the  following  equation: 

rT<.AT  -  rTrM=  kxm  (61) 


where  r™  is  the  radius  of  turn  for  <J>  ,  x  is  the  projected  miss 
distance  using  <J>r,  rTc  is  corrected  radius  of  turn,  and  k  is  a 
proportionality  constant  which  will  be  tuned  for  best  performance. 
In  other  words,  the  arc  length  is  increased  by  an  amount  propor¬ 
tional  to  the  miss  distance.  When  Eq (61)  is  rearranged,  the 
equation  for  r,pc  is 


This  law  exhibits  the  following  desirable  characteristics. 

If  x  is  zero,  then  <b  is  an  exact  solution  and  r-  equals  r™  . 
m  r  i  c  ^  1  r 

The  change  in  turn  radius  is  proportional  to  x  ,  and  inversely 
proportional  to  AT  -  the  further  away  the  tanker  is  from  roll¬ 
out,  the  less  correction  needed. 


Eq(6l)  also  gives  some  insight  into  setting  k.  A  value 

of  k=1  would  increase  the  arc  length  by  the  amount  of  x  .  Since 

m 

both  aircraft  are  flying  at  about  the  same  speed,  if  the  tanker 
would  maintain  r,pc  until  rollout,  it  would  roll  out  at  about 
the  same  time  that  the  receiver  would  reach  the  rollout  point 
based  on  .  A  good  initial  trial  value  for  k  would  then  be  1  , 
The  ad  hoc  control  algorithm  is  then: 

1.-7.  Same  as  equations  (60a)  through  (60h). 


2 


g  tan 


§rTc 


At  this  point,  the  30  degree  maximum  bank  angle  restriction 
should  be  recalled.  For  both  the  turn  point  and  the  closed  loop 
bank  angle  control  algorithms,  the  following  conditional  state¬ 
ment  was  added  after  step  2.  (Eq(60b)): 

2a.  If  <t>r  £  30  degrees,  then  set  <{>r  =  30  degrees. 

For  the  closed  loop  control  algorithm,  if  is  greater 
than  30  degrees,  <jQ  should  also  be  set  to  30  degrees  and  the 
algorithm  stopped  at  that  point.  In  that  case,  the  tanker  is 
headed  for  an  overshoot  of  the  desired  rollout  track:  it  cannot 
turn  any  tighter,  and  a  shallower  bank  angle  will  only  make  the 
overshoot  worse.  Also,  if  y ^  is  greater  than  the  tanker  has 


55 


r« 


DD 
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already  overshot  the  receiver's  estimated  cross- track  position, 
and  the  full  30  degrees  of  bank  should  be  used. 

So  far,  a  mathematical  model  of  an  air  refueling 
rendezvous  has  been  created,  and  estimation  and  control  algorithms 
have  been  developed.  The  last  major  task  of  this  study  is  to 
do  an  initial  evaluation  of  the  filter  and  controller  designs. 


?! 

9 
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Simulation  Results 


Computer  simulations  of  a  typical  rendezvous  were  used 
for  an  initial  analysis  of  the  estimator  and  controller  designs. 
Simulation 

The  filter  and  controller  were  evaluated  using  Monte 
Carlo  simulations.  The  statistics  for  each  Monte  Carlo  run 
were  calculated  from  the  results  of  30  separate,  simulated 
rendezvous.  In  each  rendezvous,  a  different  set  of  random 
variable  realizations  was  simulated. 

The  receiver's  dynamics  and  the  TACAN  measurements  were 
modelled  exactly  as  described  in  the  summary  of  Section  III. 

The  tanker's  dynamics  were  modelled  by  Eqs(4l).  (42),  and  (45). 
In  each  case,  the  tanker  started  from  a  position  along  the  minus 
y  axis  with  an  initial  flight  path  direction  parallel  to  the 
refueling  track.  vR  and  Vrp  were  both  set  to  396  knots,  which 
would  be  typical  for  a  refueling  altitude  of  25000  ft  and 
standard  day  conditions. 

SOFE  (Ref  14)  and  SOFEPL  (Ref  15).  two  computer  programs 
developed  for  Monte  Carlo  simulations  of  Kalman  filter  design 
problems,  were  used  to  accomplish  the  runs  and  generate  the 
appropriate  outputs.  The  extended  Kalman  filter  equations 
discussed  in  Section  III  are  automatically  implemented  in  SOFE 
with  one  exception.  Rather  than  the  standard  update  equations 
given  by  equations  (21)  through  (25).  a  Carlson  update  (Ref  14: 
28-29)  is  used.  This  type  of  update  is  algebraically  equivalent 
to  the  standard  form,  but  it  has  numerical  stability  and  accurac 


advantages  when  used  on  a  finite  wordlength  computer.  Whether 
those  advantages  outweigh  the  inherent,  increased  computational 
requirements  for  this  problem,  would  be  an  area  for  future  study. 

To  ensure  a  standard  means  of  comparison,  the  same 
control  algorithms  were  used  for  all  filter  evaluations.  The 
turn  point  algorithm  given  by  Eqs(60a)  through  (60i),  and  the 
cross- track  controller,  which  only  uses  Eq(56)  at  each  update 
time,  were  used.  Also,  the  control  algorithms  used  the  true 
state  values  rather  than  the  estimated  states  as  input.  For  the 
same  reason,  each  controller  was  evaluated  using  the  same  filter 
algorithm  to  provide  the  estimated  state  inputs. 

Filter 

The  big  unanswered  question  about  the  proposed  system 
is  how  much  information  can  be  gained  from  TACAN,  distance-only 
measurements.  Ideally,  the  filter  would  be  so  accurate  that 
radar  would  not  be  needed.  Then,  accurate  rendezvous  could  be 
conducted  even  if  radar  contact  is  lost  or  never  achieved.  The 
filter  would  be  started  as  soon  as  the  receiver  passes  the  initial 
point.  An  accurate  estimate  of  the  receiver's  cross- track  position 
would  then  be  used  by  the  tanker  to  establish  its  initial  offset. 
Unfortunately,  though,  the  initial  simulations  indicated  that 
attaining  such  a  degree  of  accuracy  may  not  be  feasible. 

A  7  state  filter  which  modelled  all  of  the  states  except 
the  short  correlation  time  component  of  TACAN  error  was  initial¬ 
ized  at  simulated  100  NM  ranges.  The  results  indicated  that, 
as  suspected,  some  of  the  states  were  relatively  unobservable. 


The  filter  predicted  error  standard  deviations,  as  measured  by 
the  square  roots  of  the  diagonal  elements  of  P,  decreased  very 
little  for  the  TACAN  bias,  altimeter  error,  and  constant  v 

J 

states.  Between  the  start  and  finish  of  the  rendezvous,  the 
predicted  standard  deviations  only  decreased  about  1 0%  for  b, 

for  vyC*  and  virtually  not  at  all  for  ag.  Also,  those  slight 
improvements  for  b  and  v^c  did  not  occur  until  about  the  last 
60  seconds  of  the  rendezvous. 

For  an  initial  100  NM  range  and  a  value  of  6.15  NM, 

v 

both  the  7  state  filter  and  the  3  state  filter  described  in 
Section  IV  were  very  difficult  to  tune.  The  filters  produced 
biased  estimates;  and  they  greatly  overestimated  their  own 
accuracy;  the  diagonal  elements  of  P  quickly  converged  to  values 
much  less  than  the  true  error  variances.  An  accurate  cross¬ 
track  estimate  never  could  be  achieved  until  the  range  decreased 
considerably.  In  Figure  5,  the  standard  deviation  of  the  filter's 
actual  cross-track  estimate  error,  computed  for  each  sample  time 
of  a  typical  Monte  Carlo  run,  is  plotted.  Notice  that  at  t=310 
seconds,  which  corresponds  to  about  a  30  NM  range,  the  actual 
filter  error  had  a  standard  deviation  of  around  4  NM.  The  tanker 


would  be  better  off  using  radar,  with  its  assumed  la  error  of  1 
NM  at  that  range. 

When  it  appeared  that  no  progress  was  being  made  in 
designing  an  adequate  estimator  for  that  "worst  case"  situation, 
attention  was  shifted  to  a  more  restricted  problem.  It  was 
assumed  that  the  tanker  would  use  radar  to  establish  its  offset 
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and  to  obtain  a  relatively  short  range  initial  estimate  of  the 
receiver’s  position.  The  filter  would  then  take  over  and  pro¬ 
pagate  and  update  the  state  estimate  through  the  rest  of  the 
rendezvous.  The  advantage  of  even  this  restricted  case  is  that 
it  would  provide  a  state  estimate  for  turn  point  calculations 
and  bank  angle  control  during  the  turn.  The  radar  is  good  for 
initial  estimates  and  positioning,  but  the  one  onboard  the  KC-135 
cannot  lock  on  a  target.  Therefore,  it  could  not  easily  be  used 
for  periodic  updating  during  the  turn  since  all  radar  updates 
must  be  entered  manually. 

The  filtering  was  initiated  at  a  simulated  range  of  30 
NM, giving  about  a  50  second  period  before  reaching  the  turn  point. 
This  shorter  range  also  makes  the  velocity  assumptions  about 
the  receiver  more  valid.  By  30  NM,  the  receiver  should  have 
completed  any  initial  maneuvering  or  altitude  changes  that  may 
have  been  necessary,  and  its  flight  path  should  be  fairly  stable. 

The  results  of  the  3  state  filter  simulations  are  plotted 
in  Figures  6,  7,  and  8.  All  of  the  filter  parameters  were  tuned 
to  the  values  discussed  in  Section  IV.  In  each  figure,  the 
statistics  of  the  actual  filter  error,  computed  from  the  time 
histories  of  30  simulated  rendezvous,  are  depicted.  The  actual 
filter  error  is  defined  as  the  difference  between  the  appropriate 
truth  state  value  and  the  filter  estimated  state  value.  The 
center,  solid  line  represents  the  average  error;  the  two  broken 
lines  are  the  mean  error  plus  and  minus  the  standard  deviations; 
and  the  two  dashed  lines  are  the  filter  predicted  standard 
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deviations,  calculated  as  the  average  square  root  value  of  the 
appropriate  diagonal  element  of  P. 

Each  of  the  plots  demonstrate  some  desirable  filter 
characteristics.  The  mean  error  stays  near  zero.  Also,  the 
filter  predicted  standard  deviations  fairly  closely  bracket  the 
actual  mean  plus  standard  deviation  plots.  In  other  words,  the 
filter  performs  its  updates  based  on  an  internal  model  which 
closely  approximates  its  actual  performance. 

The  plots  also  indicate,  however,  that  the  accuracy 
does  not  improve  significantly  from  the  initial  uncertainties. 

For  the  x  position  and  velocity  components,  no  noticeable 
improvement  is  made  until  about  t  =  120  sec,  which  is  when  the 
tanker  is  approximately  halfway  through  the  turn.  Meanwhile, 
the  y  estimate  just  exhibits  a  modest,  steady  improvement. 

Some  of  the  filter  parameters  were  varied,  but  no 
noticeable  improvement  was  achieved.  For  example,  shorter  sample 
periods  were  used.  The  results  of  Figures  6  to  8  were  based  on 
a  simulated  time  between  measurements  of  10  sec.  For  Figures 
9,  10,  and  11,  a  shorter  period  of  2  sec  was  used.  The  actual 
error  statistics  stayed  about  the  same. 

Another  interesting  result  was  found  when  experimenting 
with  the  simulated  TAGAN  errors.  The  TACAN  bias  error  standard 
deviation  of  0.1  is  large  compared  to  the  short  correlation 
time  component's  constant  variance  of  0.016  NM.  However,  if  the 
bias  error  is  removed  entirely  from  both  the  truth  and  filter 
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Figure  6.  Actual  and  Predicted  x  Estimate  Error  Statistics 


ALONG-TRfiCK  VELOCITY  ESTIMATE  ERROR 


Figure  7.  Actual  and  Predicted  vxe  Estimate  Error  Statistics 


Figure  8.  Actual  and  Predicted  y  Estimate  Error  Statistics 
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models,  the  filter  performance  is  still  about  the  same.  The 
cross-track  error  statistics  for  this  case  are  plotted  in  Figure 
12.  The  plots  are  almost  identical  to  those  in  Figure  8.  This 
result  indicates  that  not  much  would  be  gained  by  being  able  to 
estimate  b  or  calibrate  it  out  of  the  TACAN. 

All  of  the  simulations  seemed  to  indicate  that  the  results 
plotted  in  Figures  6  to  8  represent  a  practical  limit  for  the 
extended  Kalman  filter.  The  driving  constraint  on  filter 
performance  does  not  seem  to  be  the  number  of  filter  states,  the 
sample  period,  or  the  magnitude  of  the  TACAN  bias  error.  The 
rendezvous  geometry  itself,  would  intuitively  seem  to  be  the 
primary  constraint.  Because  the  aircraft  are  approaching  each 
other  nearly  head-on,  not  much  cross-track  information  is  avail¬ 
able  until  late  in  the  rendezvous.  Investigation  of  different 
rendezvous  geometries  might  be  a  suitable  area  for  future  study. 

The  improvement  in  estimation  accuracy  over  the  initial 
uncertainties  is  not  significant.  But  the  extended  Kalman  filter 
does  seem  to  be  a  relatively  simple,  stable  method  of  generating 
a  receiver  state  estimate  up  through  completion  of  the  rendezvous. 
With  the  present  method,  the  tanker  has  no  useful  means  of  estima¬ 
ting  the  receiver’s  state  once  the  turn  is  started. 

Controller 

The  most  important  question  of  this  study  still  remains 
to  be  answered:  what  improvement  in  rendezvous  accuracy  might 
be  achieved  with  the  proposed  method?  To  answer  that  question,  the 
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Figure  12.  y  Estimate  Error  Statistics-No  TACAN  Bias  Error 
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control  algorithms  of  Section  V  were  incorporated  into  the 
Monte  Carlo  simulations.  For  each  case  studied,  the  3  state 
filter  with  10  second  updates  provided  the  inputs  to  the  controller. 
Each  rendezvous  was  initialized  at  a  30  NM  range  and  was  terminated 
when  the  tanker  had  turned  to  the  refueling  heading.  The  rollout 
errors  were  then  calculated  as 

xm  =  x^tf)  -  2  NM  -  xT(tf) 

ym  -  y(tf)  -  yT(tf) 

where  x  and  y  are  the  rollout  miss  distances,  x  and  y  are  the 

true  receiver  position  coordinates,  and  t^  is  the  final,  or 

rendezvous  completion,  time.  For  each  Monte  Carlo  run,  30 

simulated  rendezvous  were  used  to  compute  the  means,  standard 

deviations,  and  root  mean  square  statistics  for  and  y  . 

To  form  a  basis  for  comparison,  the  present  rendezvous 

method  was  also  simulated.  For  the  airspeeds  and  zero  wind 

conditions  used  in  the  simulations,  the  turn  range  chart  (Ref  16: 

3-5)  gives  a  nominal  offset  of  7  NM  and  a  turn  range  of  17  NM. 

In  the  first  case  studied,  it  was  assumed  that  the  tanker 

had  maneuvered  to  the  computed  offset  distance,  based  on  y,  prior 

to  the  turn.  The  only  actual  offset  errors  were  then  due  to  the 

y  estimate  errors.  For  the  present,  open  loop  method  simulations, 

this  meant  that  the  tanker  had  established  an  estimated  7  NM 

offset,  but  the  actual  offset  had  a  standard  deviation  of  1  NM 

* 

due  to  radar  estimate  errors.  For  the  controlled  rendezvous, 
the  tanker  started  from  an  estimated  offset  of  9.6  NM  -  the 
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nominal  offset  for  a  25  degree  bank  turn.  Again,  the  actual 
offset  varied  from  this  mean  value  with  ala  error  of  1  NM. 

The  results  are  listed  in  Table  1.  Since  the  present 
method  has  no  means  of  improving  its  estimate  of  the  receiver's 
y  position,  its  la  value  for  ym  is  about  equal  to  the  initial 
radar  estimate  error.  Also  notice  the  large  bias  errors  at 
rollout.  These  would  primarilly  be  due  to  using  the  whole  mile 
range  values  from  the  turn  chart,  though  some  of  that  error  can 
be  reduced  by  interpolating  values  from  the  chart. 

The  cross-track  controller  only  used  Eq(56)  at  each  update 

time  to  compute  a  bank  angle  which  would  roll  the  tanker  out  on 

course.  The  1  a  ym  value  decreased  by  a  factor  of  4.  However, 

there  is  a  slight  increase  in  the  error  statistics  for  x  .  That 

m 

is  because  when  the  turn  is  started,  the  filter  still  has  about 
a  . 4  NM  y  error  standard  deviation.  As  a  better  estimate 
becomes  available  during  the  turn,  br  changes  and  so,  therefore. 


Table  1.  Rollout  Error  Statistics 


xm (NM )  : 


ym(NM) 


Controller 


Open  Loop 
Cross-Track 
Closed  Loop: 
k  =  2 
k  =  3 


RMS  | 

m 

a 

RMS 

1.09  1 

.35 

1.11 

1.16 

• 

V-O 

00 

.14 

1 

.23 

.27 

.84 

,0 

.43 

.  44 

00 

• 

1  -.10 

.41 

.42 

does  f.  The  change  in  the  turn  rate  from  that  predicted  at  the 

start  of  the  turn,  causes  the  rollout  time  and  hence  x  to  vary. 

Closed  loop  control  along  both  axes  was  simulated  next. 

Results  for  values  of  k  equal  to  2  and  3  are  listed.  With  only 

a  slight  increase  in  the  standard  deviation  of  y  ,  the  statistics 

for  xffl  improved  considerably.  Most  of  the  mean  error  is 

eliminated.  For  k=3,  the  root  mean  square  error  is  about  half 

of  what  it  is  for  the  open  loop  and  cross-track  controllers. 

Since  the  desired  values  for  x  and  y  are  zero,  their  RMS  values, 
_ _  m  J  m 

where  RMS  =  J m2+  a*  are  an  appropriate,  comparative  measure 
of  controller  performance. 

The  values  for  k  equal  to  2  and  3  were  selected  from  the 

results  of  a  larger  set  of  runs.  In  Figure  13,  the  RMS  miss 

distances  are  plotted  for  values  of  k  between  0  and  10.  The 

case  of  k  equal  to  zero  implies  cross-track  control  only:  the 

bank  angle  which  minimized  ym  was  computed,  and  no  adjustments 

were  made  to  reduce  x  .  As  k  was  increased  above  zero,  the  RMS 

m 

values  for  x  decreased,  while  the  y  RMS  values  increased 
m  *  m 

slightly.  However,  for  k  equal  to  9  or  more,  the  bank  angle 
corrections  were  too  great,  and  the  controller  produced  increas¬ 
ingly  large  rollout  errors.  At  k  equal  to  12,  the  RMS  values 

for  x  and  y  were  9.3  and  3.0  NM  respectively.  The  controller 
m  *  m 

did  not  seem  to  be  particularly  sensitive  to  changes  in  gain 
in  the  range  of  k  equals  2  to  6.  The  x  miss  distances  were 
improved  over  those  of  the  cross- track  controller,  and  the  y  miss 
distances  remained  small. 
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Figure  13.  Effect  of  Control  Gain  on  RMS  Miss  Distances 


A  Monte  Carlo  run  was  also  conducted  in  which  the 
controller  inputs  were  the  true  states  and  not  the  filter  estimates. 
For  a  gain  of  k  equal  to  3,  the  resulting  RMS  value  of  0.17  NM 
and  y  RMS  value  of  0.09  confirmed  that  there  were  no  basic  flaws 
in  the  closed  loop  control  concept  itself. 

One  other  case  was  simulated.  A  turn  from  other  than 
nominal  conditions  was  modelled.  Such  an  initial  error  might 
be  due  to  pilot  error,  insufficient  time  to  correct  to  the 
desired  offset,  weather  avoidance,  and  so  on.  For  the  open  loop, 
cross-track,  and  closed  loop  controllers,  the  turn  was  started 
from  an  offset  2  NM  wider  than  nominal,  and  5  seconds  earlier 
than  it  should  have  been. 

For  this  case,  the  differences  in  controller  performance 

were  more  pronounced.  The  results  are  listed  in  Table  2.  The 

open  loop  controller  had  large  mean  errors  because  it  never  tried 

to  compensate  for  the  initial  deviation.  The  cross-track  controller 

reduced  the  y  RMS  value  by  a  factor  of  10,  but  it  did  not  correct 
*  m  * 


Table  2.  Rollout  Error  Statistics  -  Turn  Initiated  From 
Other  Than  Nominal  Conditions 


Controller 

m 

x  (NM) 
m 

a 

RMS 

y 

m 

(NM) 

m 

a 

RMS 

Open  Loop 

.96 

.70 

1.19 

-1.22 

1.07 

1.62 

Cross-Track 

.93 

1.02 

1.38 

.04 

.13 

.14 

Closed  Loop  (k=3) 

-.15 

.74 

.75 

-.09 

.39 

.40 

75 


for  along-track  errors.  Finally,  the  closed  loop  controller  had 
significantly  lower  RMS  values  than  the  present  method  for  both 
x  and  y  .  That  was  because  the  closed  loop  corrections  reduced 
the  mean  rollout  errors  to  near  zero. 

Even  though  the  filter  does  not  supply  a  particularly 
accurate  estimate  of  the  receiver's  state,  the  closed  loop 
controller  and,  to  a  lesser  extent,  the  cross-track  controller, 
exhibit  some  advantages  over  the  present  open  loop  method. 

The  cross-track  RMS  errors  are  reduced  considerably.  The  closed 
loop  controller  also  moderately  reduces  the  RMS  along-track 
error.  Additionally,  these  control  algorithms  compute  rendezvous 
solutions  for  other  than  nominal  offsets  and  headings,  and  can 
correct  for  deviations  during  the  turn. 


VII  Conclusions 


0 


The  purposes  of  this  research  were  to  develop  a  digital 
estimation  and  control  method  for  conducting  air  refueling 
rendezvous,  and  to  evaluate  its  feasibility  and  possible  benefits. 
For  this  study,  the  design  was  constrained  by  present,  onboard 
KC-135  navigation  equipment,  and  the  basic  geometry  of  the 
current  method. 

A  3  state  extended  Kalman  filter,  which  only  incorporates 
air-to-air  TACAN  distance  measurements,  is  not  a  very  accurate 
estimator,  and  it  appears  that  a  radar  estimate  would  be  needed 
to  initialize  the  filter.  However,  it  does  appear  to  be  a  simple, 
stable  method  of  generating  a  receiver  state  estimate  throughout 
the  final  portion  of  the  rendezvous. 

Using  the  filter  generated  estimate,  a  single  equation 
control  law  can  significantly  reduce  cross-track  error  at 
rendezvous  completion.  A  control  law  can  also  be  implemented 
which  decreases  the  along-track  error  as  well.  Using  the  closed 
loop  control  algorithm,  the  root  mean  square,  cross-track  error 
at  rollout  was  reduced  by  up  to  75$  over  the  present  method; 
and  the  along-track  error  was  reduced  as  much  as  47$. 

The  proposed  method  has  two  main  advantages  over  the 
present  method.  It  can  compute  rendezvous  solutions  from  other 
than  nominal  conditions;  and  it  generates  a  receiver  state 
estimate  which  enables  the  tanker  to  use  closed  loop  control 
throughout  the  rendezvous. 
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VIII  Recommendations 


cr< 


The  computer  simulation  results  do  indicate  that  the 
proposed  system  has  some  potential  advantages  over  the  present 
method.  The  next  step  should  be  a  more  complete  analysis  of  the 
models  and  designs  developed  here.  Several  areas  for  further 
study  immediately  come  to  mind. 

The  truth  model  could  be  improved  with  a  stochastic 
model  for  the  tanker's  dynamics  as  well  as  the  receiver's. 

Such  a  model  would  be  better  for  evaluating  different  controller 
designs  because  the  effect  of  deviations  during  the  turn  could 
be  evaluated.  In  this  study,  only  initial  deviations  were 
simulated.  Such  a  model  could  also  provide  a  basis  for  a 
stochastic  control  design  analysis. 

There  are  several  estimation  issues  yet  to  be  considered. 

A  sensitivity  analysis  of  the  filter  to  modelling  errors  should 
be  done,  and  different  aircraft  types  and  airspeeds  should  be 
simulated.  The  effect  of  the  rendezvous  geometry  on  observability 
could  also  be  investigated. 

A  present  hardware  constraint  was  imposed  for  this  study. 
However,  given  the  limited  filter  accuracy  achieved,  the  potential 
advantages  of  incorporating  other  equipment  might  be  considered. 
For  example,  the  AN/ARN  139  (V)  TACAN  system  which  provides 
air-to-air  bearing  as  well  as  distance  information  could  be 
considered  (Ref  17). 
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The  potential  advantages  of  closed  loop  control  were 
investigated,  but  only  one  algorithm  was  developed.  Other  ad 
hoc  methods  could  be  investigated  as  well  as  stochastic  control 
methods. 

A  vital  area  which  was  barely  considered  in  this  study 
is  the  onboard  computer.  The  estimation  and  control  algorithms 
need  to  be  evaluated  partly  on  the  basis  of  computation  time, 
memory  requirements,  and  so  on.  Also,  the  algorithms  developed 
here  might  be  simplified  further. 

Very  little  work  has  been  done  in  the  past  on  the  air 
refueling  rendezvous  problem.  Therefore,  research  in  just 
about  any  area  of  the  problem  might  result  in  significant 
improvements. 
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